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DOMAIN E GENERAL 

L'invention concerne les techniques de traitement de donnees 
sismiques acquises au moyen de capteurs multicomposantes. 
5 Cette invention est en particulier applicable a ['acquisition au moyen 

de cables disposes sur le fond de la mer (techniques dites « OBC » ou 
« Ocean Bottom Cable » selon la terminologie anglo-saxonne generalement 
utilisee). 

II a recemment ete propose des geophones multicomposantes 
10 capables de travailler dans n'importe quelle position, notamment dans Ie 
fond de la mer, Cette technologie de sonde « toute-inclinaison » (« omnitilt » 
selon la terminologie anglaise) a permis de nouveaux cables simplifies (les 
cardans mecaniques ne sont plus necessaires) et permet des acquisitions 
avec une meilleure largeur de bande sismique. 
15 Toutefois, Petape d'acquisition ne permet pas de fournir Porientation 

reelle des geophones du cable, alors que cette information est 
indispensable pour permettre de traiter les donnees. : 

L'invention propose quant a elle un traitement qui est destined etre 
mis en oeuvre sur les donnees brutes et qui permet la reorientation ! et le 
20 calibrage (destine a convertir en une reponse commune de phase et 
d'amplitude les mesures des divers geophones). 

ETAT DE LA TECHNIQUE 

II a deja ete propose des techniques consistant a isoler du signal les 
25 donnees qui correspondent a la premiere arrivee sur le capteur et a 
determiner a partir de ces donnees un filtre destine a etre applique aux 
donnees brutes de fa?on a les corriger et obtenir ainsi les composantes du 
signal sur les axes attendus. 

Une proposition a cet effet a ete decrite dans Particle : 
30 « Horizontal vector infidelity correction by general linear 

transformation" - Joe Dellinger et al. - SEG - 9-14 Septembre 2001 . 

Toutefois, cette technique n'est pas necessairement optimale 
puisque le mecanisme de couplage qui intervient au niveau du geophone 



2 



n'est pas le meme pour les ondes qui correspondent a une premiere arriv6e 
sur le capteur et pour !es ondes reflechies ou converties par les horizons 
sismiques. 

5 PRESENTATION DE U1NVENTION 

^invention propose quant a elle une autre approche qui emploie la 
fenetre reelle de donnees pour reconstruire numeriquement des geophones 
orientes seion les axes vouius. 

Implicitement, cette approche cornpense les erreurs qui ne sont pas 
10 liees aux geophones eux-memes, mais qui sont dues au fait que le 
couplage entre le geophone et les ondes a enregistrer est different selon 
que Ton a a faire a un mouvement de vibrations vertical plutot qu'a un 
mouvement horizontal (ce du fait de la pesanteur). 

Dans le cas d'un cable, le couplage est en outre different selon que 
15 le mouvement de vibrations est dans le sens du cable ou transversal. 

Par ailleurs, puisque des fenetres plus profondes sont sujettes a un 
rapport S/N (signal sur bruit) plus faible, des traitements mettant en oeuvre 
des sommations de traces sont en outre utilises. 

^invention propose ainsi, selon un premier aspect, un procede de 
20 traitement de donnees sismiques acquises au moyen d'un capteur a au 
moins trois composantes geophone, caracterise en ce qu'on determine des 
estimateurs qui sont des combinaisons de ces composantes permettant 
. d'isoler les differentes donnees selon q^elles correspondent a une 
propagation avec reflexion ou avec conversion et en ce que, pour 
25 determiner une reconstruction de capteur, on determine les operateurs a 
appliquer aux differentes composantes du capteur de fagon a minimiser 
Tecart entre des donnees de reference et des donnees obtenues en 
appliquant les estimateurs a la reconstruction de capteur, les operateurs 
ainsi determines etant appliques aux donnees acquises. 
30 On precise ici que Pon entend, dans la suite du presente texte, par 

geophone tout capteur de vitesse et par hydrophone tout capteur de 
pression. 
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Des aspects preferes, mais non limitatifs du procede selon ie 
premier aspect de Pinvention sont les suivants : 

- Ie capteur incluant en outre un hydrophone, les donnees de reference 
pour reconstruire un geophone vertical sont derivees des donnees 

5 acquises par I'hydrophone ; 

- les donnees de reference pour reconstruire un geophone vertical sans 
hydrophone ou pour reconstruire des geophones horizontaux sont 
derivees de I'application des estimateurs a Tun des geophones du 
capteur; 

10 - ('orientation dans Ie plan horizontal d'une composante geophone est 
obtenue en minimisant I'estimateur de la reflexion transverse ; 

- les estimateurs sont determines en fonction d'un modele de propagation 
isotrope ou incluant I'anisotropie azimutale. 

Selon un autre aspect plus general, I'invention propose un procede 
15 de traitement de donnees sismiques acquises au moyen d'un capteur a au 
moins trois composantes geophone, caracterise en ce qu'on determine des 
estimateurs qui sont des combinaisons de ces composantes perraettant 
d'isoler les differentes donnees selon qu'elles correspondent .a une 
propagation avec reflexion ou avec conversion. Les estimateurs- ainsi 
20 determines peuvent trouver d'autres applications que celle faisant I'objet du 
procede selon Ie premier aspect de I'invention. 

DESCRIPTION DES FIGURES 

- la figure 1 est une representation schematique donnant les 
25 conventions angulaires utilisees ; 

- les figures 2 et 3 sont des organigrammes donnant les 
differentes etapes du traitement respectivement dans I'un et 
I'autre des deux exemples de mise en oeuvre decrits. 

30 



DESCRIPTION D'UN QU PLUSIEURS MOD ES DE MISE EN 
CEUVRE DE L'INVENTIQN 
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Premier exemple de mise en oeuvre : cas d 'un modele de 
5 propagation isotrope 

Sous I'hypothese d'une geologie localement 1D (une dimension) a 
proximite des recepteurs, et en supposant une propagation isotrope de la 
terre, un geophone donne, avec une orientation cp ip, mesure : 

m k =Rppcos(\|r) 8pp k + (Rpscos(e k -cp)+Rtrsvsin(9 k -<p))sin(\|/)8ps k 
Avec : 

k : index pour le point de tir (allant de 1 a N) 
6 k : Azimut du point de tir par rapport a I'axe des abscisses 
R PP : reflectivite PP 

6pp : Correction dynamique PP (« normal moveout » ou. NMO) 
Rps : reflectivite radiale isotrope PS 
Rtrsv: reflectivite transversale isotrope PS 
5 Ps : Correction dynamique PS (« normal moveout » ou NMO) 
Ce modele permet revaluation des parametres de reflectivite a partir 
de I'ensemble de traces tr k par un traitement simple de comparaison des 
moindres carres (en ignorant pour commencer) dans le domaine de 
Fourier, menant aux equations suivantes: 
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20 



25 



N wc((p) ws(q>) 
wc(cp) Sc2((p) Scs((p) 
ws (J) Scs(<p) Ss2(<p) 



Rps 
^Rtrsv^ 



Svx 



Shcx(9) 
Shsx(cp) 



J 



Quantites scalaires : 

30 Sc((p)=Z k cos(9 k -q>) Ss(cp)=2 k sin(9 k -(p) 

Sc2(cp)=2 k cos 2 (9 k -cp) Ss2(cp)=Z k sin 2 (0 k -9) Scs(<p)=S k cos(0 k -cp) sin(9 k -<p) 
N= Sc2+ Ss2 
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Quantites ondelettes : 

wc(<p)= E k cos(6 k -cp) 5ps k Sppk" 1 ws(cp)= Z k sin(9 k -cp) 5ps k 5pp k " 1 

Quantites sommations de traces du qeophone x : 
Svx = Z k x k 8pp k " 1 

Shcx(cp)= S k cos(9 k -cp) x k 6ps k " 1 Shsx(cp)= Z k sin(e k -9) x k 5ps k " 1 
La solution de ce systeme lineaire donne : 

Aiso Rppcos(v(/) = W Svx + (Scs ws - Ss2 wc) Shcx + (Scs wc - Sc2 ws) Shsx 
AisoRps sin(\|/) = Kic cos(cp)+ Kis sin(cp) 
Aiso Rtrsv sin(\|/) = Kis cos(q>)- Kic sin(cp) 

Avec: 

W = Sc2 Ss2-Scs 2 

Aiso = W N + wc (Scs ws - Ss2 wc) + ws (Scs wc - Sc2 ws) 

Kic = (Scs ws - Ss2 wc )Svx + (N Ss2 - ws ws )Shcx + (- N Scs + ws wc )Shsx 

Kis = (Scs wc - Sc2 ws)Svx + (» N Scs + wc ws )shcx + (N Sc2-wc wc)shsx .': 



Cette modelisation permet des evaluations en prenant en compte les 
proprietes suivantes : 
a. Rpp ne depend pas de <p, 
20 b. |Rps| 2 +|Rtrsv| 2 ne depend bien entendu pas non plus de cp, 

c. Aiso est dans la pratique rapidement stationnaire dans le temps et peut 
etre ignore pour le processus de calibrage/orientation, car commun a 
tous les geophones d'un meme recepteur. 

25 Evaluations de tirs denses 

La plupart des acquisitions OBC sont tirees en utilisant une grille 
dense et reguliere de source, qui permet beaucoup de simplification : 
Sc = Ss = 0, wc = ws=0 (symetrie des sources par rapport aux recepteurs) 
30 Scs =0 Sc2 = Ss2 = N/2 (distribution de source isotrope) 
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Ensuite, la solution exacte peut etre obtenue par I'approximation : 
N Rpp cos(\j/) = Svx 
N Rps sin(\|/) = 2 Shcx(cp) 
5 N Rtrsv sin(v|/) = 2 Shsx(cp) 

Cette approximation mene a des calculs tres simples, n'impliquant 
pas d'ondeiettes, et peut etre mise en application immediatement. 

10 Orientation des geophones 

Puisque le Rtrsv n'existe pas physiquement, la minimisation de 
Penergie de Rtrsv mene a une equation trigonometrique qui donne 
Porientation reelle cp ge0 (+k n): 

15 

i 

tg(2cpJ = 2(£Kic, Kis,)/(XKic t 2 -£Kis t 2 ) 
t t t 

((Emax-Emin)/ (Emax+Emin)) 1/2 donne un controle de la qualite de la 
reorientation. 

20 D'autre part, si I'on veut trouver I'orientation d'apres les arrivees 

premieres, il est possible de corriger lesdites arrivees premieres pour les 
mettre a un meme temps d'arrivee, puis de simplifier kic et kis en 
remplacant les ondelettes wc et ws par les scalaires sc et ss, en 
considerant que les ondes enregistrees horizontalement sont en fait la 

25 projection de I'onde P radiale, presente sur tous les geophones du fait de 
son obliquite. 

Calibrage composite vertical de qeophone : 

30 Avec les geophones g-i ,g 2 .gs, nous construisons un geophone 

composite vertical v, v = opi*gi + op 2 *g2 + op 3 *g3 (ou comportant plus de 
termes similaires dans ie cas ou des geophones supplementaires sont 
presents dans le recepteur) ou opi , op 2 , op 3 sont les filtres de longueur 
finie et op u *g u represente la convolution du geophone g u avec le filtre op u . 
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tels que : 

E1 = |XH - XV| 2 =|Kic(v)l 2 + |Kis(v)| 2 

5 L'energie de la difference entre XH (I'hydrophone apres application 

du fant6me geophone, ou hydrophone cross-ghost) et XV (le geophone 
composite vertical apres application du fantome hydrophone, ou geophone 
cross-ghost), (voir par exemple a cet effet la demande de brevet de la 
demanderesse FR 2 743 896). 

10 

E2 = |Rps(v)| 2 + |Rtrsv(v)| 2 

Energie horizontale du composite vertical, 

Alors E = AE1 + (1- A)E2 est une forme quadratique sur les 
15 coefficients des filtres et peut etre reduite au minimum, donnant ainsi un 
systeme lineaire a resoudre. (A est un parametre d'accord, 0<=A<=1, qui 
favorise soit un ajustement plus fort a I'hydrophone de reference soit une 
minimisation plus forte de l'energie de cisaillement). 

20 Dans le cas de donnees terrestres, c'est a dire s'il n'y a pas 

d'hydrophone disponible, on peut choisir un des geophones comme 
reference et remplacer I'hydrophone par Rpp(g re f). 

Calibration dans une direction arbitraire horizontale 

25 

Avec g = op^gi + op 2 *g2 + op 3 *g 3 et q> g une direction arbitraire, 
il est defini : 

E1 = lRps(g, q> g ) - Rps(ref, <p re f)| 2 , 

comme etant l'energie de la difference entre revaluation du Rps du 
30 geophone composite arbitraire et revaluation du Rps d'un geophone de 
reference (en general le geophone oriente dans la direction du cable). 

E2 = |Rpp(g)| 2 +IRtrsv(g, cp g )| 2 (l'energie non radiale) 
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E = AE1 + (1 - A)E2 permet la derivation d'un geophone horizontal 
composite dans la direction voulue, ayant la meme reponse en frequence 
que le geophone de reference, et avec une contamination PP minimale. 
5 En considerant les cas q> ref =0 et qw =7t/2 on peut generer par simple 
combinaison trigonometrique les projections radiales et transverses. 

Deuxieme exemple de mise en oeuvre : cas d' une modelisation de 
propagation azimutale anisotropique 

to 

Soit a la direction de I'axe de propagation rapide naturelle. La modelisation 
de la mesure du geophone devient (en utilisant I'un ou I'autre de Rps1 et 
Rps2 les deux images selon les directions naturelles, ou Rps et 8Rps 
definis par Rps1= Rps+SRps, Rps2= Rps-8Rps): 

15 

m k = Rpp cos(h/) 5pp k + (Rps1 cos(<p - a)cos(9 k - a) + Rps2 sin(cp - a)sin(9 k - a))sin(v|/)8ps k 
ra k = Rppcos(v) 6pp k + (Rpscos(6 k - q>)+ 5Rpscos(6 k + <p - 2a))sin(M;)8ps k 

donnant les equations normales 



20 




25 La solution de ce systeme lineaire donne: 



Aiso Rpp cos(v)/) = Inchange 

Aaniso Rps sin(v) = (Kac cos(2 a - cp)+ Kas sin(2 a - <p)) sin(2 (a - (p)) 
Aaniso Rtrsv sin(v|/) = (-Kac cos((p)- Kas sin((p)) sin(2 (a - <p)) 
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Avec: 

Aaniso = sin 2 (2 (a - (p)) Aiso 

Kac = (Sc2ws- Scs wc) S vx + (n Scs - wc ws )Shcx - (n Sc2 - wc wc )shsx 
Kas = (-Ss2wc + Scs ws) Svx + (n Ss2 - ws ws )shcx - (n Scs - ws wc )shsx 

5 

Calibration verticale de geophone composite: 

10 La methode isotrope reste applicable avec le changement 
E2=|Kac(g)| 2 + |Kas(g)| 2 

Calibration arbitraire horizon tale ou Kgnes/colonnes de qeophones 
15 composites 

L'observation de 8Rps sur le champ de donnees permet de diagnosticjuer 
la presence (ou non) d'une anisotropie azimutale significative, (la quaritite 
sin 2 (2(a~9)) 5Rps ne necessite pas la connaissance de a pour son calcul). 

20 

La methode isotrope reste applicable avec les changements 
E2=|8Rps(v)| 2 et E=X(E1+E2)+(1-X) E3. 

Des lors que a n'est g6n§ralement pas connu, un balayage sur une gamme 
25 de 7i/2 est mise en oeuvre, en utilisant la valeur de a qui minimiseEmini/EO. 
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REVEN PI CATIONS 

Procede de traitement de donn6es sismiques acquises au 
moyen d'un capteur a au moins trois composantes geophone, 
caracterise en ce qu'on determine des estimateurs qui sont des 
combinaisons de ces composantes permettant d'isoler les 
differentes donnees selon qu'elles correspondent a une 
propagation avec reflexion ou avec conversion et en ce que, 
pour determiner une reconstruction de capteur, on determine les 
operateurs a appliquer aux differentes composantes du capteur 
de facon a minimiser I'ecart entre des donnees de reference et 
des donnees obtenues en appliquant les estimateurs a la 
reconstruction de capteur, les operateurs ainsi determines etant 
appliques aux donnees acquises. 

Procede selon la revendication 1, dans lequel, le capteur 
incluant en outre un hydrophone, les donnees de reference pour 
reconstruire un geophone vertical sont derivees des donnees 
acquises par I'hydrophone. 

Proced6 selon la revendication 1, dans lequel les donnees de 
reference pour reconstruire un geophone vertical sans 
hydrophone ou pour reconstruire des geophones horizontaux 
sont derivees de ('application des estimateurs a I'un des 
geophones du capteur. 

Procede selon la revendication 1, caracterise en ce que 
I'orientation dans le plan horizontal d'une composante geophone 
est obtenue en minimisant I'estimateur de la reflexion 
transverse. 



11 



5. Precede selon Tune des revendications precedentes, caracterise 
en ce que les estimateurs sont determines en fonction d'un 
modele de propagation isotrope ou incluant I'anisotropie 
azimutale. 

6. Procede de traitement de donnees sismiques acquises au 
moyen d'un capteur a au moins trois composantes geophone, 
caracterise en ce qu'on determine des estimateurs qui sont des 
combinaisons de ces composantes permettant d'isoler les 
differentes donnees selon qu'elles correspondent a une 
propagation avec reflexion ou avec conversion. 
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Orientation de gref conriue ? 





minimisation de Rtrsv(gref)) 2 par j 
rapport a <pref 


■ = l^psfe, <Pg) - Rps(gref, <Pref)H 
E2 = :|Rpp(g.)f +|Rtrsv(g, q>g)f 





\ 

Mihimlsation de E = A£1+<1-A,)£2 pour determiner op i , op 2 , op 5 




1 Geophone horizontal oriente q> g 

\ g=Opi*gl +Op2*g2 + Op 3 *g3 



Choix de g re f | 



Orientation de g ref connue ? 




minimisation de IRtrsvCgref)! 2 par rapport a 

<Pref 



El = |R ps (g, (pg) - Rps(gref, CPref)| 2 

iE2 = |Rp P (g)l 2 +|R*sv(g, 




I Minimisation de E = XE1 + (1 - X)E2 pour determiner opi, op^^j 



Figure 2 
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Geophone vertical g = op, *g, +op 2 *g 2 +op 3 *gj 



I 



Hydrophone connu ? 



NON 



X>UI 



XH - Stack Hydrophone* fantome^eophone : 
E1 «{flpp(g* faotome^hydropfiOTie) -rate' 



E2 = \Rp$mf * |R 




Choix de gref 
£1 - |Rpp (g) - 



Minimisation de E - JUE1 + <1 -i)E2 pour determiner op 4 , op 2 , op 3 



[ Geophone vertical g = opi*gi + Qp2*g2 + Qp3*g3 
l__ 3Z - 




XH = Stack Hydrophone * fant6me_geophone I 

ie)-XHl 2 J 



El = | Rpp (g* fantdmejhydrophone 




Choix de g re f 
G = Rp P (g ref ) 
El = 1 Rpp(g) -G I 2 
1 ~ 



^ HRps(g)! 2 +|Rtrsv(g)ri 
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\ Minimisation de E = XE1 .+ (1 - X)E2 pour determiner opjTopTopJ 



Figure 3 
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